Quantum electric dipole glass and frustrated magnetism near a critical point in 

Li2ZrCu04 * 
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We report a new peculiar effect of the interaction between a sublattice of frustrated quantum 
spin-1/2 chains and a sublattice of pseudospin-1/2 centers (quantum electric dipoles) uniquely co- 
existing in the complex oxide 7-Li2ZrCu04 (= Li2CuZr04). ^Li nuclear magnetic-, Cu^+ electron 
spin resonance and a complex dielectric constant data reveal that the sublattice of Li"'"-derived 
electric dipoles orders glass like at Tg ~ 70 K yielding a spin site nonequivalency in the Cu02 
chains. We suggest that such a remarkable interplay between electrical and spin degrees of freedom 
might strongly influence the properties of the spiral spin state in Li2ZrCu04 that is close to a 
quantum ferromagnetic critical point. In particular that strong quantum fluctuations and/or the 
glassy behavior of electric dipoles might renormalize the exchange integrals affecting this way the 
pitch angle of the spiral as well as be responsible for the missing multiferroicity present in other 
helicoidal magnets. 

PACS numbers: 75.10.Pq, 76.30.-v, 76.60.-k, 77.22.-d, 77.80.-e 



The recently discovered quantum spin chain compound 
Li2ZrCu04 Ul is a strongly frustrated quantum magnet 
located on the spin spiral side of the magnetic phase di- 
agram in the very vicinity of a quantum critical point 
(QCP) to ferromagnetism. Surprisingly, contrarily to 
other related materials where ferroelectricity occurs si- 
multaneously with a spiral magnetic order yielding a mul- 
tifcrroic behavior 0;[3;|3], here, below the magnetic or- 
dering at T/v Pi 6K no multiferroicity has been found 
hitherto This is puzzling since the usually sup- 

porting features such as the low symmetry of buckled 
Cu02 chains allowing Dzyaloshinski-Moriya interactions 
and/or special positions for impurity spins are present 
in Li2ZrCu04. In addition, the unusually large magni- 
tude of the estimated main (nearest neighbor) exchange 
integral Ji ^ — 300K [l[ differs markedly from those of 
other related edge-shared Cu02 chain compounds such 
as Li2Cu02 (Ji « -230 K) @ and Ca2Y2Cu50io as well 
as from L(S)DA-|-[/ calculations for Li2ZrCu04itself 
The microscopic reason for these puzzling deviations are 
completely unclear at present. We conjecture that under 
the special conditions of a QCP the effect of additional in- 
teractions or degrees of freedom might be crucial for var- 
ious physical properties including the puzzle mentioned 
above. 

Here we report such a new peculiar feature absent 
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to the best of our knowledge in all other known chain 
cuprates which should certainly affect the quantum mag- 
netism. It is specifically related to the interaction of the 
quantum S = 1/2 spins in the Cu02 chains with the 
quantum electrical dipoles derived from tunneling Li+ 
ions. Our ^Li nuclear magnetic- (NMR), the Cu^+ elec- 
tron spin resonance (ESR) and the dielectric constant 
data show that a coupling between the active electrical 
and the magnetic degrees of freedom occurs already in the 
paramagnetic regime far above T^r. We argue that this 
peculiar effect is due to the interaction of interpenetrat- 
ing magnetic and electrical sublatticcs formed specifically 
in Li2ZrCu04. 

The orthorhombic crystal structure of 7-Li2ZrCu04 (= 
Li2CuZr04) d comprises chains formed by edge-shared 
Cu04 plaquettes running along the c-axis (Fig.[T]). Ow- 
ing to an almost 90° Cu-O-Cu bonding geometry the 
quantum 5=1/2 spins of Cu^+ ions are coupled along 
the chain by the nearest neighbor (NN) ferromagnetic 
and the next-NN (NNN) antiferromagnetic exchange in- 
teraction causing spin frustration [J]. The Cu02 chains 
in Li2ZrCu04 form planes similarly to, e.g., LiCu202. 
However, the intcrplane distance of 4.7 A markedly ex- 
ceeds that of other quasi-two dimensional (2D) cuprates. 
At variance with many other Li-containing cuprates the 
Li ions in Li2ZrCu04 occupy two different types of po- 
sitions: 46 (Lin) with a 100% occupancy and 8/ (Lij) 
with an occupancy of 50% (Fig.[T]). An unusually large 
Debye- Waller factor suggests a splitting of the Lii po- 
sition [§]. This can be modeled by a double- well po- 
tential where the Lii ion can hop over the energy bar- 
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FIG. 1: (Color online). Fragment of slightly idealized 7- 
Li2ZrCu04 structure with the supertransferred Cu-O-^Lii.n 
bonds. Arrows show schematically the oxygen ion displace- 
ments in the real structure. 

rier f/b between two contiguous 8/ positions. In the 
most pronounced quantum situation at low temperatures 
/cbT ^ Uh the Lii split site can be approximated by a 
two-level system (TLS) described in terms of tunnelling 
Li"^ ions which is equivalent to pseudospin s=l/2 centers 
with a quantum reorienting electric dipole directed along 
the a-axis: da=\e\dsz- Here d is the extent of the Lii split 
position along the a-axis. Hence, the Lij subsystem can 
be described as a 3D pseudo-tetragonal pseudospin- 1/2 
lattice with Lij planes sandwiched between the adjacent 
planes of Cu02 chains (Fig.[T|). Thus Li2ZrCu04 may 
provide a unique model system to study a synergetic in- 
terplay of two interacting low-dimensional quantum sub- 
systems: the magnetic cupratc chain planes and the Lij- 
derived electric dipole planes forming a natural laminar 
composite structure. 

The effective TLS-Hamiltonian for the Lij subsystem 
using the pseudo-spin notation can be written as follows: 

Hli = ^ ifjSizSjz + hn'^Six + ^ itjSixSjx , (1) 
i<j i i<j 

where the first term describes the NN Lij -Lij inter- 
action, the second one does the Lij tunnelling between 
contiguous 8^ positions with the frequency fl, while the 
third one docs the correlated exchange tunnelling for 
NN and more distant Lij sites. The pseudospin- 1/2 
Hamiltonian ([1]) describes also the well known transver- 
sal field Ising model frequently used in the theory of 
(anti)ferroelectrics and quantum glasses. The phase or- 
dering depends essentially on the relation between Ml 
and effective NN coupling energy /J|„, however, it is an 
inter-plane frustration and/or the coupling to the spin 
system that governs the Lij dipole ordering. Turning to 
a magnetoelectric coupling between the Cu^+ planes and 
the Lii pseudospin- 1/2 sublattice one should first note 
that the Lij ion hopping between two contiguous 81 po- 
sitions modulates the 2p orbital occupation at nearest 



oxygen ions and their displacement. The latter in turn 
modulates locally the Cu^^ crystal field, the in-plane 
and inter-plane t"*" transfer integrals within a single-band 
Hubbard-type model which describes the electron (hole) 
transfer between the CUO4 plaquettes of the magnetic 
subsystem. (Note that the magnetic Cu^+ sublattice 
can also produce an effective electric field which acts on 
the Lii dipoles.) These correlated electronic models can 
be mapped afterwards on a Heisenberg model for the 
Cu spins hence resulting in a modulation of the corre- 
sponding exchange integrals. Focussing on one of the 
Lii pseudospin sites (Fig. [1]), both for the Cu^+ crys- 
tal field Vcf and the in-plane Cu^+-Cu^+ charge trans- 
fer we deal with an anticorrelation effect leading to a 
local nonequivalence of upper and lower CUO2 chains: 
AK/(1) = -AK/(4) cx s,; AK/(2) = -AK/(3) oc s,. 
The interaction term of the two subsystems reads 

AhI = At'l {a\d2 - ald4)sz + h.c. , (2) 

where dl and dj are electronic annihilation and creation 
operators and denotes the NN-transfcr integral in the 
chain direction. The modulation of the inter-planc Cu^+- 
Cu^+ electron transfer reads as follows 

AHt^.^At^{ ald/)s.^ + h.c. . (3) 

i=l,2:j"=3,4 

Mapping the electronic part on a Heisenberg model, 
we arrive at the famous Kugel-Khomskii model [8| which 
describes the interaction of spins and pseudo-spins usu- 
ally derived from orbital degrees of freedom. It has been 
widely used to model mostly on-site interactions of spins 
and orbitals and only in the mean-field approximation. 
To the best our knowledge a case of an inter-site interac- 
tion with a glassy state of one of the components was not 
described so far. However, based on our work to be re- 
ported below, we argue that it might become relevant for 
a broad class of systems with multiple degrees of freedom 
including orbital physics, too. 

To verify the above conjectures we have carried out 
^Li NMR, Cu^+ ESR, and dielectric measurements on 
oriented polycrystals of Li2CuZr04. Owing to a small 
anisotropy of the g-factor (see below) it was possible to 
align powder particles mixed with epoxy resin in a strong 
magnetic field. After its hardening a sample with a mag- 
netic " easy" axis parallel to the a-axis has been obtained, 
as confirmed by the x-ray diffraction. 

We start with the ^Li NMR measurements that proved 
to be very useful before to study both the magnetic or- 
dering in cuprates [1, [13, and the Li ion mobility 
(see, e.g., Ref. [III). The ^Li (/=3/2) NMR spectra of 
Li2ZrCu04 samples were measured by a Tecmag pulse 
spectrometer in two orientations: H||a and H±a by 
sweeping the magnetic field at a frequency a;Ar=38MHz. 
The signal was obtained by integrating the spin-echo en- 
velope. The longitudinal T-^^ and transversal Tj"^ relax- 
ation rates were measured at the peak of the signal using 
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the stimulated echo sequence and 71/2 — -k sequence, re- 
spectively. The quadrupole splitting, typical for ''Li nu- 
clei in cuprates 1^, 12 1 and estimated to be of the order 
of 0.05 MHz, is unresolved in the spectrum which shape 
can be described by a single Gaussian line. In the para- 
magnetic state at high temperatures (T > 150K Tn) 
the ''Li NMR response of Li2ZrCu04 (Fig.[2]), seemingly 
presents a single line. However, a careful analysis of the 
lineshape and rates shows that we deal with an un- 
resolved superposition of two lines with a full width of 
0.01 T, which may be ascribed to two Li species, i.e. to 
the immobile Lin and the mobile Lij ions, respectively. 
For H±a lowering the temperature below T'^lOO K yields 
a separation of the NMR spectrum in two well resolved 
lines with a strong and unusual T-dependent frequency 
redshift and inhomogeneous broadening of the low field 
(left) line, (Fig. [2]), that can be associated with the NMR 
response of the mobile Lij ions. For H||a the signals from 
different Li sites merge and only at T > T* ~100K one 
can separate two contributions due to the narrowing of 
one of them, apparently of Lij. A characteristic tem- 
perature T* ~100K can be associated with the onset of 
the quenching of the Lii hopping between two equiva- 
lent positions on the NMR timescale, i.e. with the offset 
of the motional narrowing. It means that at T > T* 
the Lii pseudospin system reveals most likely a classical 
high-temperature paraelectric behavior. One should note 
that the T-dependence of the ''Li NMR linewidth for im- 
mobile Lin fons (high-field right line) shows up a rather 
conventional low-T rise due to a critical increase of spin 
fluctuations in the vicinity of T/v. Since the distances 
Lin-Cu and Lii-Cu are similar, one can assume roughly 
the same spin fluctuation contribution to both NMR sig- 
nals. Thus one can single out the additional contribution 
to the inhomogeneous broadening of the left line (dashed 
line in Fig.[2l3) whose T-dcpendcncc turns out to be typ- 
ical for systems with mobile Li ions (see, e.g., Ref. |11|). 
The temperature of the ceasing of the motional narrow- 
ing of the NMR line is usually correlated with the onset of 
the critical glassy freezing, which suggests T* as an upper 
boundary for the glass ordering temperature Tg ^ T* . 

The low-T behavior of Tj"^ (Fig.[2li) with characteris- 
tic values of ~ 10^ s~^ and a pronounced peak near Tm 
for both ''Li species and both flcld directions is typical 
for spin ordering cuprates [l2|. The spin-lattice relax- 
ation (SLR) rates T^^ for both Li species reveal a very 
similar T-behavior with nearly the same constant value 
of ^ 750sec~^ down to low T ^ lOK where both start 
to increase with a peak at T ~ 6 K for the ''Lij line and a 
divergence for the ''Lin hne in the case of H_La. A similar 
behavior is observed for the H||a orientation where con- 
tributions from ''Liii and ''Lii nuclei are hardly resolved. 

A typical mobility induced SLR rate for ''Li NMR 
in a wide number of nonmagnetic solids is due to the 
quadrupole mechanism which is usually small as com- 
pared to a strong SLR due to a magnetic mechanism. In 
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FIG. 2: (Color online), (a) - Selected ''Li NMR spectra of 
the oriented Li2ZrCu04 sample at ujn ~ 38 MHz for H_La; 
(b) - T-dependence of the NMR linewidth for low- and high- 
field NMR lines. Dashed curve denotes an inhomogeneous 
broadening due to a glass-like ordering of Lii ions. Inset shows 
the behavior of resonance fields; (c) and (d) - T-dependences 
of the spin-lattice T^^ and spin-spin relaxation rates Tj"^, 
respectively. In (b) - (d) O denote data for H||a, ■ and • 
are data for low- and high-field lines for H±a, respectively. 



Li2ZrCu04 where strong magnetic fluctuations take place 
in the anisotropic lattice the decrease of relaxation rates 
below 100 K reflects the ceasing of the motion that was 
averaging local magnetic fields created by fluctuations 
of Cu magnetic moments. The anisotropic character of 
this motion becomes apparent in the different relaxation 
behavior in different field orientations. One can see in 
Fig.[2t,d that Ti for H || a increases concomitantly with 
T2 for H ± a, which are both caused by magnetic fluc- 
tuations perpendicular to the a-direction. To summarize 
this part, all ''Li NMR data cleary indicate a freezing of 
the Lii paraelectric sublattice below T* < 100 K. 

A direct evidence of a glass-like structural ordering of 
this sublattice can be provided by measurements of the 
dielectric constant e = e' + ie" that have been performed 
with pressed pellets of Li2ZrCu04 at a number of fre- 
quencies /. The behavior of the real part e' above 300 K 
(Fig. [3]) evidences most likely a contribution of Lii to ionic 
conductivity. However, the most strong change of e' oc- 
curs in the T-range 50- 150 K. A significant reduction of 
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FIG. 3: (Color online). Temperature dependence of the real 
(e') and imaginary (e") parts of the dielectric constant for a 
pressed sample of 7-Li2CuZr04 at different frequencies /. T- 
dependence of the remanent polarization P,. is shown by open 
circles. Solid line is a guide for the eye. 



e' clearly indicates a 'randomly antiferro' freezing of the 
paraelectric subsystem in Li2ZrCu04 [isj . Remarkably, 
in this T-range the electric field dependence of the po- 
larization P{E) exhibits a hysteresis (Fig. 2^) with the 
values of the remanent polarization (Fig.[3]) compara- 
ble with those of multiferroic cuprates [3|, |j| . 

The frequency dependence of the step in e'{T) con- 
comitant with the shift of the peak temperature T^"max 
of the imaginary part e"{T) to lower T with decreas- 
ing the frequency / (Fig. [31 inset, and Fig. 5^) reveals a 
typical for a glass-like transition critical dynamic slow- 
ing down towards the glass transition at Tg. The tem- 
perature Te"max Can bc associated with the frequency / 
dependent freezing temperature below which the longest 
relaxation time of the system gets larger than the char- 
acteristic observation time 1// and the system enters a 
noncquilibrium regime In the three dimensional 

(3D) case where a finite glass ordering temperature Tg 
can be expected, / and Tp/z^m, can be connected by the 
power-law relation / = fat^'^ [15|. Here /o is the fiuctu- 
ation frequency, t = {T^"max/Tg — 1) is the reduced tem- 
perature, and z and v are critical exponents. The / vs. 
Te"max dependence can be well fitted by this expression, 
yielding the glass ordering temperature Tg ~ 70 K and 
zv ~ 5.5, Fig.|4)D. The value of Tg relates well with the 
characteristic NMR 'freezing' temperature T*. Notably, 
the obtained product of the critical exponents zv ~ 5.5 
is very close to the result zi^ = 6 ± 1 of a reported Monte 
Carlo simulation for the 3D Ising glass [16|. This sup- 
ports the appropriateness of the Ising (pseudo-spin-1/2) 
model adopted above for the description of the electric 
Lii subsystem in Li2ZrCu04. 

How can the magnetic Cu^+ subsystem respond to this 
effect? First of all one may anticipate a change of the 
crystal field and of the g-factors for Cu^+ ions which can 
be detected by the ESR technique. We have carried out 
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FIG. 4: (Color online), (a) - Hysteresis loop of the electric 
field dependence of the polarization P(E) at T = 55 K. Inset; 
enlarged view around the origin with the remanent polariza- 
tion Pr indicated by the arrow. The temperature dependence 
of Pr is shown in Fig.|3l (b) - Full circles - relation between 
the measurement frequency / and the reduced temperature 
t = (T'e"max/?g — 1). Solid line - fit to the power law re- 
lation / = /ot^^with the values of the fluctuation frequency 
/o = 1.9 MHz, the glass ordering temperature Tg ~ 67.4 K 
and the critical exponent zv — 5.45 (see the text). 



high-field ESR (HF-ESR) experiments at sub-THz fre- 
quencies in a broad temperature ran ge w ith a home-made 
spectrometer (for details, sec Ref. [17|)- In the high-T 
regime above Tg ~ 70 K the ESR spectrum comprising a 
single Lorentzian absorption line has been found (Fig.[5l 
inset). The resonance field corresponds to g-factors of 
g|| = 2.19 and g±=2.02 for H|ja and H_La, respectively, 
typical for a Cu^+ ion in a distorted square planar ligand 
coordination [3] ■ A remarkable evolution of the ESR 
spectrum has been observed upon cooling the sample be- 
low Tg where the spectrum begins to split into two lines. 
The shape of the spectrum can be perfectly fitted with 
two Lorentzians of comparable intensity, (Fig. [5]), inset. 
Representative T-dependences of the resonance fields at 
a frequency of 352 GHz obtained from the fit are shown. 
Fig. El The two lines in the HF-ESR spectrum whose 
separation increases with decreasing T can be straightfor- 
wardly assigned to the occurrence of two nonequivalent 
Cu sites with slightly different g-factors. Since the g- 
factor is very sensitive to the symmetry and the strength 
of the ligand electrical field potential Vcf , the splitting of 
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FIG. 5: (Color online). Inset: HF-ESR spectra at 352 GHz for 
H _L a at 100 K and 15 K, black solid lines. Dash lines (red 
online) are fits indistinguishable from the experiment. De- 
composition of the spectrum at 15 K in two lines is shown by 
dash (blue online) and dash-dot (green online) curves. Main 
panel: T-dependence of the resonance field of the peak(s). 



the spectrum evidences that the mobile Lij ions freeze in 
the lattice on a particular pattern yielding two distinct 
local crystal field potentials Vcf±AVcf at the Cu^+ sites. 
The shift of both lines to higher fields indicates the on- 
set of low-frequency spin correlations which develop in 
a low-D spin system far above T/v. Owing to a much 
shorter timcscalc of the ESR they are detected at higher 
T compared to the NMR relaxation data (Fig. [2j:,d). 

Finally, we argue that the frozen 'glassy' Lij-derived 
quantum dipoles may produce a random electrical field 
potential that might disturb the oriented 'multiferroic' 
dipoles induced in some way by the magnetic spiral with 
the aid of special impurity spins or by the Dzyaloshinski- 
Moriya interaction. As a result no multiferroicity will 
be observed below the magnetic ordering at Tm- Fur- 
thermore, the absence of any long-range order of the 
Lij-dipoles comprising in Li2ZrCu04 a regular (although 
frustrated) electrical dipole sublattice is noteworthy, be- 
cause a glassy ground state was observed so far only for 
essentially nonstoichiometric random systems [isj . The 
only feasible source for an effective randomness at Lij 
sites in our high-quality samples should then be due to 
short-range spiral-like correlated spins. Note that even 
a perfectly ordered but intrinsically incommensurate spi- 
ral state is expected to introduce a quasi-random poten- 
tial at the commensurate Lii positions. Thus, one con- 
jectures that owing to a sufficiently strong coupling of 
spins and pseudospins short-range spin correlations seen 
by ESR already at T ^ Tg ^ Tjv may introduce disorder 
also in the electrical Lij-dipole sublattice. An additional 
decrease of e' below 30 K concomitant with a rapid 
decrease of Pr (Fig. [31) may indicate a crossover from a 
classical regime dominating the response of the glassy 



ordered Lij sublattice at higher temperatures to a quan- 
tum tunneling regime which interferes with the develop- 
ment of quantum magnetism in Cu02 chains. In this 
context it would be of considerable interest to suppress 
or to modulate the electrical dipole disorder, e.g., by ex- 
ternal electric fields or pressure, and to investigate how 
the paraspiral state in between Tg and Tjv as well as 
the spiral state below T/v will response. In particular, 
in the region Tn < T < Tg the renormalization of the 
chain exchange integrals due to the spin-pseudo-spin in- 
teraction might take place. As a consequence, it would 
require different fits for the magnetic susceptibility x{T) 
[H above and below Tg, being thus possibly helpful to 
resolve or to contribute to the solution of the above men- 
tioned Ji-puzzle for Li2ZrCu04. This way, our observa- 
tion of a pseudo-spin (electric dipole) - spin-1/2 coupling 
paves a way to new experimental and theoretical studies 
of Li2ZrCu04 to get more insight into the fundamental 
features of multiferroicity, magnetic ordering, as well as 
the nature of the critical point itself. 

In conclusion, by measuring '^Li NMR, Cu^+ ESR and 
dielectric constants of the quantum spin-1/2 incommen- 
surate chain cuprate Li2ZrCu04 we have observed a new 
peculiar effect of the interaction between spins- 1/2 as- 
sociated with Cu^+ ions and pseudospins- 1/2 (Ising-like 
electrical dipoles) due to tunneling Li"*" ions. From the 
NMR and dielectric measurements we conclude that at 
Tg~ 70 K a glass-like structural ordering of the electri- 
cal sublattice takes place. Its remarkable impact on the 
spin system is evidenced by the development of spin site 
nonequivalency at T < Tg revealed by ESR. These re- 
sults put forward Li2ZrCu04 as new type of the intrinsic 
single-phase magnetoelectric composites where quantum 
magnetism of Cu02 chains meet the electrically active 
regular sublattice of Li"*" ions. As such it emerges as 
a unique model system to study the influence of addi- 
tional interactions and degrees of freedom on frustrated 
magnetism near a critical point where the magnetic spin 
subsystem is especially soft. One of the feasible conse- 
quences of such influence could be, e.g., the suppression 
of multiferroicity in the title compound. 

Support of the DFG (grants 436 RUS 113/936/0-1, 
DR269/3-1 and KL1824/3-1) and the RFBR (grants 
08-02-91952-NNIO-a, 07-02-01 184-a, 06-02-17242, 07-02- 
96047, and 08-02-00633), is gratefully acknowledged. YA 
acknowledges support of the EU Programme Alban. We 
thank A. A. Zvyagin and N.M. Plakida for discussions. 



^ Electronic address: 'v.kataev@ifw-dresden.de' 
[1] S.-L. Drechsler et ai, Phys. Rev. Lett. 98, 077202 (2007) 
[2] M. Fiebig, J. Phys. D: Appl. Phys. 38, R123 (2005) 
[3] Y. Naito it et al., J. Phys. Soc. Jpn. 76, 023708 (2007) 
[4] S. Park et al, Phys. Rev. Lett. 98, 057601 (2007) 
[5] Y. Tarui et ai, J. Phys. Soc. Jpn. 77, 043703 (2008) 



6 



[6] W. E. A. Lorenz et al., in preparation S.-L. Drechsler et 

al., J. Phys: Conf. Ser. to be published (2009) 
[7] M. Schmitt et al, Phys. Rev. B 80, in press (2009) 
[8] K. I. Kugel and D. I. Khomskii, Sov. Phys. Usp. 25, 231 
(1982) 

[9] C. Dussarrat et al., J. Sol. St. Chem., 166, 311 (2002) 
[10] A. A. Gippius et al, Phys. Rev. B 70, 020406 (2004) 
[11] D. Brinkmann et al., Phys. Rev. B 26, 4810 (1982) 
[12] C. Kegler et al. Phys. Rev. B 73, 104418 (2006) 
[13] A. Loidl, Annu. Rev. Chem. 40, 29 (1989); U. T. Hochh 
et al. Adv. Phys. 39, 405 (1990) 



[14] K. Binder and A. P. Young, Rev. Mod. Phys. 58, 801 
(1986) 

[15] P. C. Hohenberg and B. I. Halperin, Rev. Mod. Phys. 49, 
435 (1977) 

[16] A. T. Ogielski and I. Morgenstern, Phys. Rev. Lett 54, 
928 (1985) 

[17] C. Golze et al, Phys. Rev. B 73, 224403 (2006) 
[18] J. R. Pilbrow, Transition Ion Electron Paramagnetic Res- 
onance, Clarendon press, Oxford, 1990 



